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MODELING AND ENGINEERING SIMULATION OF 
NUCLEAR THERMAL ROCKET SYSTEMS 

. Modular Thermal Fluid Solver with Neutronic Feedback 

. Main Component Modules: 

Pipes, Valves, Mixer 
Nozzle Skirt 
Pump, Turbine 
Reflector, Reactor Core 

. Hydrogen (Para- and Dissociated) Properly Package 

10 S T < 10,000 K 
.1 i Ps 160 bar 

• Models Developed for NTVR, NERVA and XNR 2000 
. CFD and Heat Transfer Models for Main NTR Components 


A detailed program for modeling of full system nuclear rocket engines is developed. At 
present time, the model features the expander cycle. Axial power distribution in the 
reactor core is calculated using 2- and 3-D neutronics computer codes. A complete 
hydrogen property model is developed and implemented. Three nuclear rocket systems are 
analyzed. These systems are: a 75.000 Ibf NERVA class engine, a 25,000 Ibf cermet fueled 
engine and INSPl's nuclear thermal vapor rocket. 
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NUCLEAR THERMAL ROCKET SIMULATION SYSTEM 


unit* gyit*»s. 

KRSYSOl 

KRSTS02 

KXSY503 


A jjl£*gS rROriATlES 



T0X1IH1 UQPUL1S 

SIKFLiriCD TRIAMOOYNAMIC 
SINGLE STAGE KZAW-UKE ANALYSIS 
TWO STAG! MRAN-LINE ANALYSIS 


co>x rtQDULEs 

KRATT * WHITNEY ROCKET 
NUCLEAR VAROR TKIMCAL ROC KIT 
NEAVA DERIVATIVE ROCKET 


KQl lLZ MODUL E 

ISIHTRORIC WITH IDS* COEFFICIENTS 


The main program links ail the component modules and iterates to arrive ar the user 
specified thrust chamber pressure and lemperntme and thrust level. Reactor power and 
propellant flow rate are among outputs of the simulation program. Fuel elements in the 
core module are prismatic with variable flow area ratio. Each module divides the relative 
component into N segments. 
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Axial temperature distribution of NVTR fuel surface and propellant in an average power 
rod. Reactor power is adjusted to achieve the thrust chamber temperature and pressure 
of 2750 K and 750 psi. respectively. 



Normalised axial power distribution in C C composite fuel matrix NTVR, calculated by 
DOT-2 S„ code. The axial power shape factor is an input for the simulation code. 
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Specific Impulse vs Chamber Pressure 
INSPI-NTVR @ 75000lbf Thrust 



Parametric study of thrust chamlx*r pressure and temperature impact on Isp of N WR. At 
higher pressures Isp is less sensitive to thrust chamber temperature. 
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Turbine Pressure Ratio vs Chamber Pressure 
INSPI-NTVR @ 75000lbf Thrust 



Turbine pressure ratio is sensitive to both thrust chamber pressure and temperature. For 
thrust chamber pressure of 1200 psi and temperature of 3000 K, the turbine pressure ratio 
of 1.26 is well within the range of available technology. 
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NERVA Coro Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 



Axial temperature profiles for NF.RVA- 75,000 Ibf engine are presented. The maximum fuel 
temperature is 3490 K at .7 m from the core entrance. 
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P&W XNR2000 Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=25000lbf 



Normalized Length 


Axial temperature distribution in XNH 2000 core is presented. XNR 2000 features a two 
path folded flow core fueled with (HMIVI*. The maximum fuel temperature is 3000 K at 
about 85% from the entrance to the inner core region. 
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INSPI-NTVR Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 



INSPI-NTVR Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 
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Mach Number 


INSPI-NTVR Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 



INSPI-NTVR Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 
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Pump Pressure Rise vs Chamber Pressure 
INSPI-NTVR @ 75000 1 bf Thrust 



Chamber Pressure (psi) 


Specific Impulse vs Chamber Pressure 
INSPI-NTVR @ 75000lbf Thrust 



Chamber Pressure (psi) 


c-a 
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Core Exit Mach # vs Chamber Pressure 
INPSI-NTVR <§> 75000lbf Thrust 



Chamber Pressure (psl) 


Turbine Pressure Ratio vs Chamber Pressure 
INSPI-NTVR @ 750001b! Thrust 



Chamber Pressure (psi) 
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Turbine Blade Speed vs Chamber Pressure 
INSPI-NTVR @ 75000lbf Thrust 



Turbine Blade Diameter vs Chamber Pressure 
INSPI-NTVR <g> 75000lbf Thrust 
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Axial Power Shape Factor 


P&W XNR2000 Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=25000lbf 



P&W XNR2000 Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=25000lbf 
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P&W XNR2000 Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=250001bf 



P&W XNR2000 Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=25000lbf 



NP-TIM-92 


649 


NTP: Systems Modeling 





Axial Power Shape Factor Temperature (K) 


NERVA Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 



NERVA Core Axial Flow Profile 
Tc = 2750K Pc = 750psi F=75000lbf 
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Pump Pressure Rise vs Chamber Pressure 
NERVA @ 75000lbf Thrust 



Chamber Pressure (psi) 



The hydrogen property generator utilizes two interpolation techniques and a least-square 
curve fitting routine with a pentad spline function which links least-square fitted pieces 
together. The property generator package is incorporated into the N PR simulation code 
and also into a system of CFD-HT codes. 
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Cp Versus Temperature for 
Para- and Dissociated Hydrogen 



Ar higher temperatures, the hear capacity data displays smooth lwhavior. The sharp 
increase in C p value at temperatures al>ove 2000 K is due to hydrogen dissociation. 
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Cp Versus Temperature for 
Para- and Dissociated Hydrogen 



Heat capacity of hydrogen near the critical point shows large gradient and oscillatory 
behavior. At p - 2.35 MPa the property package indicates a sharp peak for C p . 
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Thermal Conductivity Versus Temperature 
for Para- and Dissociated Hydrogen 



The hydrogen property package is a combination of two subpackages covering the 
temperature ranges 10 - 3000 K and 3000 - 10,000 K, respectively. The large change of 
gradients in hydrogen viscosity at 3000 K indicates a non physical flaw in the model. 


Cp Versus Temperature for 
Para- and Dissociated Hydrogen 



Temperature (K) 
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Cp (j/kfl-K) 


Cp Versus Temperature for 
Para- and Dissociated Hydrogen 



Thermal Conductivity Versus Temperature 
for Para- and Dissociated Hydrogen 


K 

(W/m-K) 



Temperature (K) 
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K 

(W/m-K) 


r 


Thermal Conductivity Versus Temperature 
for Para- and Dissociated Hydrogen 



Temperature (K) 


— ■ ' insti 

University of Florida 

NUCLEAR DESIGN ANALYSIS PACKAGE 

. Multigroup Cross-sections Generated by COMBINE (ENDFB-V) 

. MCNP (4.2) for Complex Geometries 

. BOLD VENTURE (3-D, Diffusion) for Power Profile and 
Reactivity Calculations 

. ANISN (1-D, S n ) for Analysis of Heterogeneous Boundaries 

. DOT rv (1, 2-D, S„) for Analysis of Reflector 

. XSDRNPM (1-D, S n ) TWODANT (2 D, S n ), NJOY, AMPX for 
Cross-comparison 
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Pressure Drop Correlation Comparison 



Nusselt Number (Nil) and Pressure Drop ( AP) 
Correlation Comparison with CFD Analysis 

(I) CFD Analysis 
Energy Equation 


0 ( 


(e + a x )u + r yi v - 


A 

dX 

Numerical Algorithm: 
ume method 


Conductive Meat Flux 


0T\ d ( . . ,,dT\ _ 

h 'dx) + W{ lt + 0 ’ )V+T ’’ U -'''dY)=*' 

MacCormack hybrid implicit •explicit, Unite vol* 



Convective Heat Flux 


S = M'A. - n) 

f pC p uTdA 
T* / pCpvHA 

A 

Convective Heat Transfer Cocffient 


Nusselt Number 


A, 


*c[gL, 

r w - 


Nu s 


A C D 

A' c 
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n DIFFUSION APPROXIMATION 


- -+■*> - -^?r VT ■ - X ' VT 




Using ihc perfect gas law. 


*r 


jjq^ T 4 

3«pbP 


WHERE 

*« 

K 

P 

T 


Rossebnd Mean Opacity 

Stefan^ Boltzmann Constant 

Photon Collision Cross Section per Molecule 

Boltimann’s Constant 

Gas Pressure 

Gas Temperature 


□ APPROXIMATION BY USING 1-D EQUATION OF RAMATLYE TRANSFER 

r 

i(r) * j‘a I (r'l exp ( -a (r-r 0 ) dr' 
q » -q # - <\ * X (i 4 (r> ■ » (r)) 


WHERE \soltdwin 

i 4 (r) : Radiarion Intensity in the Positive Direction( From Gas to Boundary) 
i (r) : Radiation Intensity in the Negative Direction( From Boundary to Gas) 

I(r) : Source Function (»aT*/n) 


imaginary 
alane in 


q' gas 



q 


Nnssclt number & Prnndtl number 


He = 
Pr = 


pul) 


riWi) 

MnKyn) 

K t m 


(111) Pressure Drop 
Compressible Flow 




,lL±h 

2 

P ] 


Incompressible Flow 




/ = 0 0014 f -/?c~ 
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(II) Nustclt Number Correlations 


(1) Colburn Equation 

Nu = 0.023/?r°*M 

(2) Dittus-Boclier Equation 

A'ii =5 0.023/?e° ^Pr 0 3 

(3) Sicdcr-Talc Equation 

( V 0 M 

(4) rttukov Equation 

»- 5 ir© 

.V = I 07 + 12.7(/vi - 0 ( 2 ) 

/ = 00014 f '- nr 0 ™ 

8 

(5) Kamu»n-BoeIter Maritnflli Equation 


*. = - 


for r 


0.833 + 5/n(5Pr + 1) + 2. 


JL 


«"(«<=#)) 


/= 0.0014 + ™ 

8 


Axial Distance Correction 


Nv{x) 


= /V» ^1.957^1 -f 



/Vu(j) — Nu 


(■ 


+ 





Figure 6.2 Velocity distribution for a fully developed turbulent 
flow in lube. (Re«1.6 E+4) 
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Heat Flux. MW/M**2 



2/0 


Figure 6.17 Nussolt number va. axial position for a developing 
isothermal pipe flow at a Reynolds number of 53000. 
(Nu) u , is the Husselt number evaluated at z/D-20. 



Z/D 


Figure 6.15 lleat transfer rates obtained by Navier-Stokes solver 
for various boundary cell slxe. A 60x00 grid ia used. 
<T„-4000 K, T w - 1800 K ( Fi*«l atm, and F^,,»0.5 atm| 
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Z/D 

Figure 6.16 Heat transfer rates obtained by Navicc-Stokes solver 
for various boundary cell sire. A 60x00 grid is used. 
(T^-4000 K , T„“X0OO K, Pi.-lO atm, and P^-9.5 atm) 



Z/D 


Figure 6.19 Comparative result between diffusion approximation and 
1-0 Integral approximation for varying the gas opacity 
due to different flow conditions. 
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